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On the basis of a microscopic model of self-consistent field, the thermodynamics of the many-partide Fermi 
system at finite temperatures with account of three-body interactions is built and the quasiparticle equations 
of motion are obtained. It is shown that the delta-like three-body interaction gives no contribution into the self- 
consistent field, and the description of three-body forces requires their nonlocality to be taken into account. 
The spatially uniform system is considered in detail, and on the basis of the developed microscopic approach 
general formulas are derived for the fermion's effective mass and the system's equation of state with account 
of contribution from three-body forces. The effective mass and pressure are numerically calculated for the po¬ 
tential of "semi-transparent sphere" type at zero temperature. Expansions of the effective mass and pressure 
in powers of density are obtained. It is shown that, with account of only pair forces, the interaction of repul¬ 
sive character reduces the quasiparticle effective mass relative to the mass of a free particle, and the attractive 
interaction raises the effective mass. The question of thermodynamic stability of the Fermi system is consid¬ 
ered and the three-body repulsive interaction is shown to extend the region of stability of the system with the 
interparticle pair attraction. The quasiparticle energy spectrum is calculated with account of three-body forces. 
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1. Introduction 

The self-consistent field model is an effective approach for describing systems of a large number 
of particles, even in the case when the interaction between particles cannot be considered to be weak 
and the density of a system to be low. This model is applicable both for systems with a finite number of 
particles, such as atomic nuclei |T}{3) or electronic shells of atoms and molecules @1IS), as well as for many- 
particle systems, when describing them by methods of statistical physics. The method of self-consistent 
field is especially useful in studying spatially non-uniform systems and phase transitions (§][7). The self- 
consistent field model serves as an efficient main approximation in constructing a perturbation theory, 
including such a theory for spatially non-uniform systems [Q and many-particle systems with broken 
symmetries urn 

Systems of a large number of particles, obeying Fermi statistics, within a phenomenological approach 
at large enough densities are described in the language of the Fermi liquid theory fl2l 1T3I . which origi¬ 
nally was developed for normal limitless systems. Afterwards, the Fermi liquid theory was generalized 
both for systems of finite dimensions |[l4j [15] and for many-particle Fermi systems with broken phase 
symmetry possessing superfluid and superconducting properties fl6l . The Fermi-liquid approach is, in 
essence, a phenomenological variant of the self-consistent field theory 11711181 . 

In contrast to the Fermi-liquid approach, in which the interaction of quasiparticles is described phe¬ 
nomenologically by means of interaction amplitudes, the account for the interaction between real parti¬ 
cles is laid in the basis of the developed microscopic approach. Usually the interaction between particles 
is described by means of pair potentials, on the assumption that the presence of other particles does 
not influence the interaction between the two selected particles. Meanwhile, for particles possessing an 
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internal structure, the interaction between a pair of particles is changed due to the presence of a third 
particle, that can be taken into account by introducing potentials depending on the coordinates of three 
particles. Such a representation follows from the consideration of exchange and multipole interactions 
of more than two particles in different orders of the perturbation theory 11911218 . Essential is the fact that 
contribution from three-body forces not only gives quantitative corrections to characteristics of a system 
calculated with account of only pair interactions, but can be necessary for qualitative understanding of 
some effects. 

Accounting for three-body interactions is also important in the theory of nuclear forces (3] [15. ■22), 
because they model the dependence of the nucleon-nucleon interaction potential on density. In particu¬ 
lar, in the case of the interaction proposed by Skyrme [23] [24], three-body interactions are described by 
a simple delta-like potential. The role of three-body forces within the framework of the quantum field 
approach is discussed in reference [25l . 

The effects of three-body interactions should also manifest themselves in the interaction of structure¬ 
less particles located in a polarizable medium, for example electrons in a lattice. However, this issue is 
completely unexplored so far. 

In this paper, self-consistent field equations are obtained for normal (non-superfluid) Fermi systems 
at finite temperatures with account of three-body forces within the approach developed earlier for Fermi 
systems with pair interactions l9l lT0llT7llT8l . In general, the mean (or self-consistent) field approxima¬ 
tion can be formulated in different ways and, accordingly, different variants of this approximation are 
possible. In most cases this approximation is introduced at the level of equations of motion. In the sta¬ 
tistical description of many-particle systems it is convenient to introduce the self-consistent field model, 
most consistently and in most general form, on the level of Hamiltonian rather than on the level of equa¬ 
tions of motion. This permits in a natural way not only to derive the self-consistent equations of motion, 
but also to build the thermodynamics of a many-particle system already within the scope of a particular 
model. 

If the self-consistent field model is formulated in a way such that all thermodynamic relations hold 
exactly in it, as it takes place for the ideal quantum gases (26), then the structure of the model and all 
its parameters will be determined uniquely. The method of constructing a model within the Hamiltonian 
approach consists in splitting the exact initial Hamiltonian into a sum of an approximating Hamiltonian 
and a residual correlation Hamiltonian. The approximating (or self-consistent) Hamiltonian is chosen as 
a quadratic operator form of most general kind, that contains indeterminate fields at the initial stage. 
Such a Hamiltonian, describing a many-particle system in the self-consistent field approximation, is re¬ 
duced to a diagonal form and the fields entering it are determined from the variational principle. Within 
this description, the concept of quasiparticles naturally arises in the microscopic approach. The effects 
associated with the interaction of quasiparticles are accounted for by the correlation Hamiltonian and 
can be calculated using the perturbation theory. It is worth noting that this method is so general that it 
can be also used for a consistent description of normal and superfluid Bose systems [27] El), in particular 
phonons (29), as well as relativistic quantized fields with broken symmetries I30II31I . 

From the derived self-consistent equations accounting for both pair and three-body interactions, it fol¬ 
lows that the delta-like three-body interaction gives no contribution into the self-consistent field. There¬ 
fore, a description of three-body interactions requires accounting for their nonlocality. Thermodynamic 
properties of the spatially uniform system are studied more in detail in the paper, and general formulae 
are derived for the quasiparticle effective mass and the system’s equation of state at finite temperatures 
with account of both pair and three-body interactions. It is shown that, with account of only pair forces, 
the repulsive interaction reduces the quasiparticle effective mass relative to the mass of a free particle, 
and the attractive interaction raises it. 

The effective mass of fermions and the pressure of system are numerically calculated for the “semi¬ 
transparent sphere” potential at zero temperature, and expansions of these quantities in powers of den¬ 
sity with account of three-body interactions are obtained. The influence of three-body forces on the stabil¬ 
ity of the many-particle Fermi system is considered and it is shown that accounting for three-body forces 
of repulsive character extends the region of stability and can lead to stabilization of the system with pair 
attraction between particles. The quasiparticle energy spectrum is calculated with account of three-body 
forces. 
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2. Hamiltonian of the Fermi system with account of three-body 
interactions 

Potential energy of a system of N particles possessing an internal structure can be represented as a 
sum of pair, three-body, etc. interactions 

C/(ri,r 2 ,... > rjv) = X^r i ,r J 0+ £ U{r i} Vj,v k ) +..., (2 1} 

i<j i<j<k 

where [7 (r j , ry) = [/(ry.r,-), U{ri,rj,r k ) is a symmetric function in all permutations of its coordinates. In 
the second quantization representation, the Hamiltonian of the many-particle system with account of 
pair and tree-body interactions has the form 


H=T+V 2 + V 3 . (2.2) 

Here, 

T = fdqdq’'¥ + {q)H 0 (q,q')'¥(q') (2.3) 

is the kinetic energy and the energy in external field Uoiq), and 

H 2 

H 0 (q, q') = -—A8{q-q') + U 0 (q)8(q- q'). (2.4) 

2m 

The energies of pair V 2 and three-body V 3 interactions can be written in the form 

= dqdq’'¥ + {q)'¥ + {q')mq,q')'i'(q')'i'{q), (2.5) 

V 3 =^fdqdq’dq"'¥ + {q)'i> + (q')y + (q")U(q,q',q")'t’(q")'i'{q')'i>{q). (2.6) 

In the formulae Ml and below, the symbol q = (r, cr) designates the space coordinate r and the 
spin projection a. We assume s = 1/2. The field operators obey the known anticommutation relations 
|26| . In what follows, we consider the case when all the potentials depend only on the space coordinates: 
U{q,q') = U (r, r'), U{q,q', q ") — U (r, r', r") . Besides that, we assume that external potential does not de¬ 
pend on the spin projection, so that Uo(q) - t / 0 (r) . In order not to take into account the condition of 
conservation of the total number of particles in all computations, we assume that the considered sys¬ 
tem can exchange particles with a thermostat. To account for it, a term with the chemical potential q is 
introduced into Hamiltonian, so that Hoiq, q') in |2.3| is replaced by one-particle Hamiltonian 

H[q,q') = H 0 (q,q')~ q8{q-q'). (2.7) 

In what follows, when talking about Hamiltonian, we imply that it includes the term with the chemical 
potential. 


3. The self-consistent field model with account of three-body 
interactions 

In order to proceed to the description of a many-particle Fermi system within the self-consistent field 
model, let us represent the total Hamiltonian HI as a sum of two terms, i.e., the self-consistent Hamilto¬ 
nian Hq and the correlation Hamiltonian He- 


H = Hq + H c . 

The self-consistent Hamiltonian is defined by the relation 

H 0 = Jdqdq'y + {q)[H{q,q') + W{q, q')}V{q') + E 0 , 


(3.1) 


(3.2) 
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which is quadratic in the field operators of creation and annihilation. Equation J3.2) includes the self- 
consistent potential W{q, q'), representing the mean field acting on a single particle, as well as the non¬ 
operator term Eq, both still indeterminate. Owing to hermiticity of the Hamiltonian, the next property 
holds W{q, q') - W* {q 1 , q). Note that taking account of the non-operator term in > 3.2 1 is essential for a 
consistent description of the thermodynamics of a system within the considered approach. The correla¬ 
tion Hamiltonian accounts for all the effects, not accounted for in the self-consistent field model: 


He = - 


^Jdqdq , 'l' + {q)'¥ + (q')U(r,r')'¥(q')'¥(q) 

^ J dqdq'dq" W + (q)'¥ + (q')'¥ + (q")U(r, r' ,r") y ¥ (q")'? (q')¥ (q) 
J dqdq’ V + {q)W(q } q’)^{q’) - E 0 . 


(3.3) 


Obviously, the total Hamiltonian H has not changed in consequence of the performed decomposition. 
Transition to the self-consistent field model consists in describing the system using the approximate 
quadratic Hamiltonian instead of the exact Hamiltonian o Entering l |3.2} , so far indeterminate, 
quantities W(q, q’) and Eq should be chosen in an optimal manner, as it will be done below. The effects 
conditioned by the correlation Hamiltonian can be accounted for by the perturbation theory flllQHTll. In 
this paper we confine ourselves to consideration of only the main approximation. Note also that in con¬ 
structing the self-consistent field theory, we do not take into account the effects connected with breaking 
the phase symmetry and those leading to the properties of superfluidity and superconductivity (9] JT0l|. 

Since Hamiltonian ( |3.2) is quadratic in the field operators, it can be represented in the form of Hamil¬ 
tonian of an ideal gas of quasiparticles. The field operators and the operators of creation and annihilation 
of free particles at, a ; are connected by the relations 

J J 


'P(fl) - I> f(<7) a jt ¥ + (<7) - 2>' 0) *(<7) a], 

j i 


(3.4) 


where the functions cf>^\q) are solutions of the Schrodinger equation for free particles 


/ 


d q’H 0 {q, q^^iq’) = e^ 0) (p ( °\q). 


(3.S) 


Here, j = (v, cr), where v is a full set of quantum numbers describing the state of a particle except the 
spin projection a. To represent Hamiltonian J3.2) in the form analogous to that for free particles, let us 


(3.6) 


introduce the quasiparticle operators yt, jj connected with the field operators > 3.4 1 by the relations 

- EWr;. ^ + (<7) = L<P*jWJr 
j j 

where the functions <pj{q) are now solutions of the self-consistent equation 


J dq'[H{q,q') + W(q,q')](pj{q') = ej(pj{q). 


(3.7) 


Note that for cpdq ) functions, the same as for (p i0 \q ), the conditions of orthonormality and completeness 
hold J 

f dq(p*{q)(l>j'(q) = 8jj’, E^y Wf ^ = 8 ^~ </)• (3.8) 

J i 

As a result, the self-consistent Hamiltonian {12} acquires the form of the Hamiltonian of an ideal gas of 
quasiparticles 

Ho = E e ;7y7y + £o, (39) 

i 

where Ej means the quasiparticle energy and Eq means the energy of the background on which the 
quasiparticles exist. In the self-consistent field model, the reference of the energy cannot be arbitrary 
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but, as will be shown below, should be chosen in a way such that the thermodynamic relations are satis¬ 
fied. From l |3.4|) and (376) it follows that the quasiparticle operators are explicitly expressed through the 
operators of free particles 



(3.10) 


In a spatially uniform non-magnetic state, the operators of particles and quasiparticles coincide, which 
corresponds to the known reasoning in the Fermi liquid theory regarding the invariance of classification 
of states during adiabatic “switching-on” of the interaction | T21[13 ;|. 

4. Derivation of the self-consistent potential 

Let us define the statistical operator 


Po = exp/1(Q - Hq), 


(4.1) 


where /3 = \IT is the inverse temperature, the constant Q = -rin[Spe“^ ff °] is determined from the nor¬ 
mality condition Sppo = 1 and, as will be seen below, has the meaning of the thermodynamic potential of 
the system in the self-consistent field model. The average of an arbitrary operator A in this approximation 
is defined by the relation 


(A) = Sp(poH). 


(4.2) 


Since Hamiltonian (379) is quadratic, then for the averages with the statistical operator (4.1) the Bloch-de 
Dominicis (Wick) theorem (32) holds. 

The ground state energy in the self-consistent field model is determined from the requirement of 


equality of the averages calculated with the statistical operator (4.1) for the exact GED and for the ap¬ 
proximating GEU Hamiltonians: 


(H) = (H 0 ). 


(4.3) 


Hence, we have 



- Jdqdq'W(q,q')('¥ + (q)'V(q')). 


(4.4) 


Let us define the one-particle density matrix by the relation 


p(q.q') = (V + Wmq)) = 5>/(<7)<W)/i- 


(4.S) 


Here, the quasiparticle distribution function is defined by the expression (yt jj) = fiSjj. Based on the 
form of Hamiltonian (3.9), a straightforward calculation gives the Fermi-type distribution function 


fi = /(£() = [expire i) + 1] 1 


(4.6) 



Since the quasiparticle energy £, is the functional of f), the formula (4.6) represents a complicated non¬ 

linear equation for the distribution function, being similar to that which takes place in the Landau phe¬ 
nomenological theory of a Fermi liquid (T2). The energy Eq expressed through p(q, q [ ) has the form: 



3 p{q, q)p[q', q”)p{q", q') + p(q, q)p(q', q')p(q", ?")] 



(4.7) 
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The variation of the thermodynamic potential Q = -rin[X„(n|e P H °\n)] is equal to the averaged 
variation of Hamiltonian (3.2) : 


dn = 


T.n(n\ 

e " 

-P h °8H 0 \ 

n) 

Ln( 

n\ 

\ e -P H o\ 

\ n 

) 


(8H 0 ). 


(4.8) 


From t he r equirement that the variation of the thermodynamic potential with respect to the density 
matrix 4.5 ' vanishes 8Q/8p(q, q') = 0, we obtain the expression for the self-consistent potential which 
consists of the contributions from pair and three-body interactions W{q,q') = W^ 2 \q,q') + W^\q,q'), 
where 


W i2 \q, q') = -U{v,r')p(q, q') + 8{q- q') jdq"U (r, r")p(q", q"), 


(4.9) 


W {3 \q, q') = J dq" U(r,r',r”) [p(q, q")p{q", q’) - p(q, q’)p{q", q")} 

+ l -8{q- q’) J dq"dq’" U{ r,r", r'")[p{q", q")p(q’", q'") - p{q ", q"')p{q"', q")\. (4.10) 

The self-consistent field can be also derived by an equivalent method from the requirement that the 
variation of the functional I = (H - Hq) with respect to the density matrix vanishes, where W{q,q') 
and the parameter Eq are not varied (I1ED3- Such variation rule is conditioned by the fact that, as may be 
checked, the variations of the functional 81 with respect to 8W ( q, q') and 8Eq are mutually compensated. 
Equation s together with the derived potentials (4.9) , (4.10) , enables us to obtain the quasiparticle 
wave functions cpjiq ) and the quasiparticle energies ep. 

e i = -^jdq<l)* i (q')A(pi(q) + jdqU 0 {r)\(pi(q)\ Z + Jdqdq' W(q, q')(p*(q)cpi(q') - p. (4.11) 

Substitution of the potentials (4.9) , (4.10) into equation (3.7) leads to the integro-differential equation 

H 2 

-—A(pi(q)+ [ t/ 0 (r) — /i — £;] (pi [q) 

2m 

+ |J dq 1 U{r,r')p(q', + dq’dq" U(r,r',r")[p{q', q')p{q", q") - p{q', q")p(q", q')\ 1 (piiq) 

+ Jdq'dq" [7(r,r',r") [p{q, q")p{q", q') - p{q, q')p(q”, q")]ipi{q') 

-J dq' U{r,r')p(q,q')(p i {q') = 0. (4.12) 

The chemical potential p is associated with the average number of particles N by the relation 

N = / drn(r), n(r) = ^p(t 7 ,^). (4.13) 

J O 

In many cases, finding the equilibrium characteristics of the researched system does not require the 
calculation of the quasiparticle wave functions, but it is sufficient to know the one-particle density matrix. 
From equations (4.5) and (4.12) , the equation for the one-particle density matrix follows 

ft 2 

— [A p{q, q') - A 'p(q, q')] - [C7 0 (r) - Uo(r')] p(q, q) 

+ J dq"[U(r,r") - U(r',r")][p(q, q")p{q", q') - p(q", q")p{q, q')] 

+ Jd^M^q[/(ry\r"V[/(i y ,r^r''0][p(^^p(^^^0-p(^^0p(^^^ , )]p(^,^) 

+ \piq,q') $dq"dq w [U{ry y y')][p{q" ,q")ptq"\q'")-ptq\ (4.14) 
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In the absence of magnetic effects p{q,q') = p{r,r')8 (T(T i, and the self-consistent potential is diagonal 
in spin indices as well W{q,q') = [W t2) ( r,r') + W t3) (r,r , )]<5 (7(7 ', where 

W (2 \ r,r') = - U (r, r')p(r, r') + 28{r- r') J dr" U (r, r") p (r", r"), (4.15) 

W l3) (r,r') = Jdr"Lf(r,r',r")p(r,r")p(r",r')-2p(r,r') Jdr" [/(r,r',r")p(r",r") 

+S[ r - r') J dr" dr"' U(r, rV") [2p (r", r")p(r"', r'") - p (rV")p(r"\r") ]. (4.16) 

In what follows we consider the system without account of magnetic effects. 


5. Thermodynamic relations 


Let us formulate the self-consistent field model at finite temperatures in a way such that all the ther¬ 
modynamic relations would hold already in this approximation. This requirement leads to a unique for¬ 
mulation of the model. Entering the definition of the statistical operator quantity O has the meaning of 
the grand thermodynamic potential in the self-consistent field model. The entropy is defined through the 
statistical operator {4.1} by the known expression 

S - -Sp(p 0 lnpo). (5.1) 


With the help of ED it is easy to verify that the usual thermodynamic definition of the grand thermo¬ 
dynamic potential holds D.-E - TS - pN, where E is the total energy of the system. This potential is a 
function of the temperature T and of the chemical potential p, as well as a functional of the one-particle 
density matrix p(q, q’) - p(q, q’; T,p), which also depends on these quantities. However, by virtue of the 
fact that the self-consistent potential was derived from the condition 8El/8p{q, q’) = 0 , when finding the 
derivatives of O with respect to temperature and chemical potential, one should account for only the 
explicit dependence of the thermodynamic potential on these quantities. As a consequence, at a fixed 
volume of the system, the usual thermodynamic relation proves to be fulfilled in the self-consistent field 
model: 

dO = -Sdr-ATdp. (5.2) 


Calculation of the entropy and the number of particles either by means of the averaging with the density 
matrix {4.1} or with the help of the thermodynamic relations S = -(dQ IdT)^, N = -( dClldp) T gives the 
same result. With the correct choice of the energy Eq, any kind of inconsistency in statistical description 
of systems within the self-consistent field model, mentioned in Hi, does not appear. 

The energy (4/7), with account of the obtained self-consistent potential {4.9}, {4.10}, acquires the form 


E 0 = i Jdqdq 1 U(r,r')[p(q,q')p(q',q)-p[q,q)p(q’,q')} 

+Jdqdq'dq"U (r,r',r ")\p(q, q)p(q', q")p{q", q')-^p(q, q)p{q', q')p{q”, q") 

-^p(< 1, q')p(q', q")p{q" , p)], (5.3) 


and the thermodynamic potential is determined by the formula 


n = £ 0 - r£ln(l + e _/3£ ']. 
i 


(5.4) 


The total energy can be calculated either by averaging the Hamiltonian operator or by the use of the 
thermodynamic potential {5.4} 


E = 


Q-p 


dO dO 

- T — 

dp dT 


(5.5) 
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and has the form 


E = J^Eifi +pN+^Jdqdq' U(r,r') [p[q, q')p{q', q) - p(q, q)p(q', q')] 

+Jdqdq'dq"U(r,r',r"){p(q,q)p(q',q")p(q",q')-^p(q,q)p(q',q')p(q",q") 

~p(<7, q')p{q', q")p(q", q )]. (5.6) 


The entropy in this model is expressed through the distribution function formally in the same way as 
in the model of an ideal gas 

S = - £ [fi In/; + (1 - ft) ln(l ( 5 . 7 ) 

i 

but, as remarked previously, the distribution function itself is derived from the complicated nonlinear 
equation (4.6) and includes the effects conditioned by both pair and three-body interactions. 

The total number of particles l|4.13|l is expressed through the distribution function by the formula 


N=ll ft. 

i 


(5.8) 


Since, owing to interaction, the energy of a single particle ( |4.11| l differs from the energy of a free particle 
and includes collective effects, then in this case a particle should be treated as a quasiparticle and the 
function f] as the distribution function of quasiparticles. According to (4.13) , l |5.8} , the number of initial 
free particles and the number of quasiparticles coincide, as it takes place in the Landau theory of a normal 
Fermi liquid (12). In the presence of pair correlations, leading to the property of superfluidity, the number 
of quasiparticles is less than the number of initial free particles GSQ39, because some fraction of particles 
participates in the formation of the pair condensate. 

Within the framework of the microscopic model of self-consistent field, formulae describe 

in general form the thermodynamics of the many-particle Fermi system at finite temperatures with ac¬ 
count of pair and three-body interactions. 


6. Interaction potentials in the self-consistent field model 


In calculations within the self-consistent field approximation, instead of realistic potentials by which 
individual particles interact in the vacuum and that usually model a strong repulsion at short distances, 
the effective potentials with a set of adjustable parameters are used (23] [24). A delta-like potential is 
often chosen as the simplest model potential of the interparticle interaction. Such a choice of the pair 
interaction potential U{ r',r") = t 2 8[r' - r") results in the following form of the self-consistent field 

W (2f (r,r')=|i5(r-r')n(r), (6.1) 

where n{ r) = 2p(r,r) is the particle number density. In this case, the contribution into the self-consistent 
field (4.15) of the first (exchange) term proves to be two times less by absolute value than the contribution 
of the second term, describing the direct interaction, and has an opposite sign. 

As can be immediately verified, when choosing three-body forces in the delta-like form 


!/(r',rV") = f 3 <5(r' - r")<5(r' - r'"), 


(6.2) 


which is used, for example, in (23), the self-consistent potential W t3) (r,r') 4.16 vanishes. Therefore, in 
order to obtain a nonvanishing contribution of three-body forces into the self-consistent field, we have to 
take into account their nonlocality, that is their finite radius of action. Physically this property is condi¬ 
tioned by the Fermi statistics, not allowing for the presence of more than two particles with opposite spins 
in the same space point. In the case of an arbitrary three-body potential of the interparticle interaction, 
its contribution into the self-consistent potential (4.16) can be represented in the from 


W (3) (r, r') = ai(r, r') - 2p (r, r') a 2 { r, r') + <5 (r - r') [2 bi (r) - b 2 (r) ], 


(6.3) 
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where 


ai(r,r') = J dr" Lf(r,r',r")p(r,r")p(r",r'), 
b,(r) = J dr'dr" U(r,r',r'')p (r', r')p (r", r"), 


a 2 ( r,r') = J dr" [/(r,r',r")p(r",r"), 

b 2 (r) = Jdr'dr" (/(r,r , ,r")p(r , ,r")p(r",r'). (6.4) 


Let us discuss the issue of choosing the potentials of interparticle interactions, which are suitable for 
the self-consistent field model. In a spatially uniform state, the interaction potentials should satisfy two 
conditions. They must not change under the replacement r —<• r + a (a is an arbitrary vector) of all the 
vectors they depend on, and must be symmetric relative to any permutation of their arguments. The pair 
potential satisfying these conditions has the form 


t/(r,r , ) = i[f/ 2 (r-r') + [/ 2 (r'-r)], 


(6.5) 


where U 2 (r) is the function of a vector argument. The three-body potential, depending on the pairs of 
vector differences and satisfying the similar conditions, is as follows: 

U( r,r',r") = - [f/ 3 (r-r',r-r") + t/ 3 (r-r",r-r') + U 3 (r' -r,r'-r") 

6 

+ U 3 (r - r", r' - r) + U 3 (r" - r, r" - r') + U 3 (r"-r',r"-r)]. (6.6) 

Here, U 3 {r,r') is a function of two vector arguments. If there is a dependence only on the distances be¬ 
tween particles, and the function in 16 . 61 is symmetric U 3 (r,r') - U 3 {r',r), then C7(r, r') = C/ 2 (|r — r'|) and 


f7(r,r',r") = ^[t/ 3 (|r-r'|,|r-r"|) + f7 3 (|r'-r|,|r'-r")| + [/ 3 (|r"-r|,|r"-r'|)]. 

In particular, the three-body potential can be chosen in the form proposed in reference 11211 : 

|r-r'| + |r-r"| 


yf|r-r'|,|r-r"|) = Uq exp - 


ro 


There is also a possibility of choosing the potential in the form of the Gauss function: 


t/(|r-r'|,|r-r"|) = 


u 0 


r 3/2 r 3 


exp 


'0 


(r-r) z + (r - r r 


(6.7) 


(6.8) 


(6.9) 


Such a choice is characteristic because in the limit r 0 — 0, the potential |6.9| turns into the potential 
of zero radius (6.2) , and its contribution into the self-consistent field vanishes. In principle, the model 
three-body potential can be chosen to depend on the three distances between three particles 


U{r, r', r") = U 3 (|r - r' |, |r ■- r"|, |r' - r" |). 


( 6 . 10 ) 


Note that the derivation from first principles of the potential of interaction of three atoms as structureless 
entities presents a complex quantum mechanical problem (33). 

The presently developed programs use different types of the model interatomic three-body potentials. 
These potentials are often of a complex form and in most cases include the angular dependence (3411 , 
essential in the formation of the crystalline phase. In our paper, which can be used, for example, for the 
description of the system of 3 He atoms in the gaseous or liquid phases, we restrict ourselves to the use of 
potentials of the simplest form, not depending on angles. 

The potentials which rapidly tend to infinity at small distances are often used for describing the in¬ 
teraction between particles. An example of such a potential is given by the Lennard-Jones potential 


U L] (r)=4 £ 



( 6 . 11 ) 


containing two parameters: the distance a and the energy e. The use of the potentials with the hard core 
leads to considerable difficulties, especially in the quantum mechanical description (35). Such potentials 
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do not have the Fourier image, and the self-consistent field becomes infinity when using them. Some¬ 
times, this fact is used as justification of non-applicability of the self-consistent field model in this or that 
case, for example for describing the liquid (36). Thus, the issue of the choice of the interparticle inter¬ 
action potentials is of principal significance for problems of statistical physics. The use of the potential 
which rapidly tends to infinity at small distances means that an atom or another compound particle re¬ 
tains its identity at arbitrary high pressures. Meanwhile, it is clear that the critical pressure should exist 
at which the atoms approach each other so closely that they should be “crushed” and loose their identity. 
Therefore, the requirement of absolute impermeability of particles at arbitrary high pressures is exces¬ 
sively rigorous and unphysical and, in our opinion, it is more reasonable to use the potentials which 
tend to a finite value at small distances. An example of such a potential is given by the known Morse 
potential Uyiir) = e{exp[-2(r- Tq)! a] -2exp[-(r- Tq)I a]\, where e is the energy parameter and rg, a are 
specific distances. Notice that, from the quantum mechanical point of view, the use of bounded potentials 
means the possibility for particles to tunnel through each other with some probability. The problem of 


calculating the integral jU(r )dr diverging for the potentials of type ' 6.11 1 is encountered, for example, 
when deriving the known Gross-Pitaevskii equation, which is presently widely used in describing atomic 
Bose-Einstein condensates (37). Here, the diverging integral is replaced by a finite value coinciding with 
the scattering length under the fulfilment of the Born approximation, which in essence means the use 
of the potential not diverging rapidly to infinity at small distances. It should also be noted that quantum 
chemical calculations lead to the potentials having a finite value of energy at zero |38ll'39l . 

The model of “semi-transparent sphere” potential is quite often used as the simplest form of the pair 
potential where the noted problems with divergences are absent: 


U 2 {r) = 


U 2m , r < r 2 , 
0 , r > r 2 . 


( 6 . 12 ) 


This potential is determined by two parameters, one of which U 2m has the energy dimension and defines 
the interaction strength, and the second one r 2 of the length dimension defines the radius of interaction. 
Such an observable quantity as the scattering length in the Born approximation (which, as should be 
noted, is far fro m bei ng always valid for realistic potentials) can be expressed through the parameters 
of the potential' 6.12 : ao = (, mlAnH 2 ) / U 2 [r)dr = mU^r^/Sh 2 . A similar model can be also used in the 


case of three-body forces. For the potential 
have 

U 3 [r, r') = 


6/7) depending on pairs of distances between particles, we 


U 3 m , 
0 , 


r < r 3 , r' < r 3 , 
else. 


(6.13) 


Here, there are also two parameters: the strength f/ 3m and the radius r 3 . For this choice, the total potential 
f/(r,r',r") = U 3m , if the distances between each pair of three particles are less than r 3 . However, if two 
distances between particles are less than r 3 , and the third one is greater than r 3 and less than 2r 3 , then 
f/(r,r',r") = (7 3m /3. In other cases, the potential vanishes. 

For the potential (6.10) , depending on three distances between particles, the model of “semi-trans¬ 
parent sphere” potential is defined by the formula 


l/ 3 (r,r',r") = | 


U 3m , r < r 3 , r' < r 3 , r" < r 3t 
0 , else. 


(6.14) 


In this case, the potential is nonzero only under the condition that the distances between each pair of 
three particles are less than r 3 . Let us consider more in detail the spatially uniform system of Fermi 
particles in the absence of an external field and with account of three-body forces. 


7. The spatially uniform system 

In the spatially uniform state, the one-particle density matrix is a function of absolute value of the 
coordinate difference p(r,r') = p(|r-r')|, and the particle number density n = 2p(0) is constant. For the 
potential (6.7) depending on pairs of distances between particles, the coefficients (6.4) which determine 
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the contribution of three-body forces into the self-consistent potential acquire the form 


fli(r-r') = ^Jdr"U 3 {\r-r'\,r”)p(r")pt\r-r‘ 
a 2 { r-r') = Jdr" f/ 3 (|r-r'|, r") + ^ fdr" U 3 (\r-r' + r"|, r"), 


■' + r" I) + ^ J< dr" U 3 (\r-r' + r"\, r") p(r") p (|r - r' + r" |), 


bi = p 2 ( 0 ) 


Jdrdr'[/ 3 (r,r'), 


h 2 = ^ Jdrdr' f7 :i ([r- r'|, r')pV) + | drdr' 


I /3 (t, r')p 2 (r'). 


(7.1) 


Expansions of the three-body potential and the one-particle density matrix in Legendre polynomials can 
be used: 


U 3 {\r-r'\,r") = £ [/ 3 ,(r,r';r")P/(cos0), 


z=o 


p(|r-r'|) = ]T p/Cr.r'lPKcose), 

z=o 


(7.2) 


where 

1 

Usi{r f r';r") = J r 2 + r' 2 -2rr’x, dx, 

-1 

1 00 

2 I _|_ j _ 2 1 _|_ ^ n 

p/(r,r') = —I p(\/r 2 + r ' 2 - 2 rr'xjP/(x)dx = J f{£ k )jiUcr)jiUcr')k 2 dk, (7.3) 
-1 0 


ji{x) is the spherical Bessel function. With these expansions taken into account we have: 


fli(r) = 

8 n r 
- mr, 


■J J 

0 

a 2 {r) = 

OO 

87tp(0) r 


0 

hi = 

16^ 2 p 2 (0) 

— 

32ji z r 2 


, , ,9 , , 4tt “ 1 

')p 0 {r,r')r' 2 dr' + — Y 
6 /=< 


OO 

00 1 r 

;? 0 2^1 I L 


U 3l (r, r'\r')p(r')pi(r, r')r' 2 dr', 


mr y )r '2 Arl+ 


OU 

/ 


U 3 q{t, r'; r')r' 2 dr', 


OO OO 


0 0 


’.4) 


In the spatially uniform case, the state of a single particle can be characterized by its momentum, and 
equation (3.7| admits solutions in the form of plane waves: 


, , , & 00 ' ikr 


(7.5) 


where q = (r,cr), j = (k,cr'). As was noted, in the absence of magnetic effects p{q,q') — S aa ip( r-r'), 
W(q, q') - 8 aa 'W( r-r'). Equation j4.11 1 gives the quasiparticle energy not depending on the spin projec¬ 
tion: 

h 2 k 2 

£k=- - p+W k , (7.6) 

2 in 


where W k - f drW(r)e lkr . The self-consistent potential can be represented as a sum of direct and ex¬ 
change terms 

W (r) -WqS (r) + Wc(r), (7.7) 
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and, if the interaction depends only on a distance between particles, then also W c (r) = W c (r). In this case, 


21.2 


£k- 


W-k 

2m 


4n 


CXJ 

/‘ 


■ H+Wo+ — I dr rWc(r)sin(kr). 
/c 


(7.8) 


Considering that both pair and three-body interactions give contribution into the self-consistent potential 
W 0 = W 0 (2) + W 0 (3) , W c (r) = W [2 \r) + w£ 3) (r), where 


W 0 l2) - nUm, 
W n (3) - 2b] - b 2 , 


W [2 \r) = -U 2 {r)p{r), 
wj?\r) - a\{r) - 2p{r)a 2 (r). 


(7.9) 


Here, C/ 2 (r) is the potential of the pair interaction, U 2 o = fU 2 (r) dr, n = 2p(0) is the particle number 
density, and quantities a\{r), a 2 {r), b\, b 2 for the potential (6.7) are defined by the formulae (7.4) . The 
one-particle density matrix has the form 


OO 

P(r) = -J dfcfcsin(fcr)/(£jt), 

o 


(7.10) 


where /(£fc) = [exp(/3e fc ) + l] _1 . 

At low temperatures in a degenerate system, its properties are determined by quasiparticles located 
near the Fermi surface. In this case, the notion of the effective mass of a quasiparticle m, can be intro¬ 
duced. We define the Fermi wave number at finite temperatures by the relation 

b 2 k 2 4 n 

- H+Wq +— / drrWc(r)sin(A:Fr) = 0. (7.11) 

2m Ay J 

o 


Then, near the Fermi surface k = Ay + Ak, so that the dispersion law of a quasiparticle can be represented 
in the form ty = {H 2 k?lm*)/S.k, where the effective mass is defined by the formula 


1 1 An r , 

- = j drr 3 W c (r)Mk ¥ r), (7.12) 

m* m J 

o 


and j i(x) = (sinx- jccosx)/x 2 is the first order spherical Bessel function. It is seen that the effective mass 
is determined by the exchange part of the self-consistent potential. It should be stressed that within the 
self-consistent field model, the formula (7.12) is exact, and it is fair for any densities at which the use 
of the self-consistent field approximation is still permissible. Within the scope of the model itself, the 
constraint on density cannot be derived. For estimation of the limiting value of density at which the self- 
consistent field model remains true, one should calculate a correction due to the correlation Hamiltonian 
0 by the perturbation theory and requires it to be small relative to the main approximation (lT)- The 
contribution of three-body interactions into the effective mass of a quasiparticle with different choices 
of the model potential will be studied in detail in a separate paper. In the approximation of the effective 
mass, the distribution function acquires the form 


fk- 


exp 


h 2 k? 

P~ — (k- kp) 


m » 


+1 exp 




2m* 


(7.13) 


+ 1 


For degenerate systems, both representations of the distribution function (7.13) are equivalent. The only 
difference between the functions (7.13) and the distribution function of an ideal gas is the dependence of 
the effective mass entering (7.13) on temperature and density. 

The effective mass can be introduced for the non-degenerate system as well, when the main contri¬ 
bution into the thermodyn amic s is given by particles with small momenta. Having used the expansion 
sinfcr « kr- (fcr) 3 /3!, from 1 7.8 1 , we find the dispersion law of quasiparticles in the non-degenerate case: 
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e k - h 2 k 2 l2m* c - p*, where the effective mass m* c and the chemical potential p» are now determined 
by the formulae: 


OO 

- - 

m* c m 3 J 
o 

OO 

/ 


r 4 W c (r), 


/ drr Wc(r). 


(7.14) 


Note that the formula ' 7.12 1 for the effective mass turns into the formula 1 7.14 1 , if jiik^r) « k^rl 3 is put 
in f TV2\ . 

The thermodynamic potential per unit volume, with account of three-body interactions depending on 
pairs of distances between particles and in the approximation of the effective mass, according to js3J, 
45.41 has the form 


Q 


2 T 


77 - ~J3 


V 


h 2 k 2 
2m „ 


+ 47T 


uu <_xj 

J (7 2 (r)p 2 (r)r 2 dr-2p 2 (0) J U 2 (r)r 2 dr 


+167T 


OO OO 

)Jp 2 (r)r 2 drJ[ 




)r ,2 dr'+ -p(0) / p 2 (r)r 2 dr / U 30 (r, r'; r')r'*dr' 


p(0) Jr 2 drf U 3 {r, r')p 2 (r' 

0 0 0 0 
DO OO OO OO 

P 3 8»pdr/ W r,rtr»d^-|/p(r,^d r / (% (r,rV(rV.(r,rtr»dr 


(7.15) 


Here, A= [2nh 2 1m t T) 1/2 is the thermal wavelength. The first term contains one of the integrals of the 
Fermi-Dirac function tabulated in the McDougall and Stoner paper (40) 

OO _ 

z s 1 dz 


® s (f) = 


4 f: 


• + 1 


(7.16) 


[F(s) is the gamma function], by which all the thermodynamical quantities of an ideal Fermi gas can be 
expressed. The first term in (7.15) gives the contribution of a gas of non-interacting quasiparticles. Its 
difference from the thermodynamical potential of an ideal Fermi gas consists only in the replacement of 
mass by the effective mass. The rest terms in (7.15) are specific for the self-consistent field model and are 
determined by pair and three-body interactions. Since the pressure is associated with the thermodynam¬ 
ical potential by the known relation p = -D./V, then the formula 1 7.15 ' determines, except for sign, also 
the pressure. The formula for the particle number density 




h 2 k 2 \ 

(i - 

y 2 m, 


(7.17) 


together with (7.15), determines the pressure as a function of density, that is the equation of state. 


8. Fermi system at zero temperature with three-body interactions in the 
model of "semi-transparent sphere" 

The system of Fermi particles at zero temperature can be investigated in detail if the three-body inter¬ 
action is chosen in the form (6.10) which leads to somewhat more simple formulae than for the potential 
(6.7). For the potential (6.10), the parameters determining the self-consistent potential (6.4) are given by 
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the formulae 


oo oo 

OO 1 n r 

ai(r) = 47 r^ ——- U 3l (r,r')p/(r,r')p(r')r' 2 dr', a 2 (r) - 47 tp( 0 ) U 30 (r, r')r' 2 dr', 

1 =0 ZI + 1J J 


OO oo 


hi = 16?r 2 p 2 (0) 


0 0 


oo oo 


t/ 3 o(r,rV 2 dr', 


h? = 167+ 


W 
0 0 


U 30 (r,r’)p 2 (r')r' 2 dr'. (8.1) 


Here, 


U 3 {r,r',\r-r'\) - ^U 3l ir,r')Pi{cosd), 
l=o 

l 

21+1 


2 ^ j n 

U 3l (r, r) = ——— / 1 / 3 ( 7 -, r', \J r 2 + r ' 2 - 2 rr’x)Pi(x) dx. 
-1 


*0 


where x 0 = (r 2 + r ' 2 - r|)/2rr'. At zero temperature fa - 0(fc F - fc) and 


Taking into account the latter relations, we find 


( 8 . 2 ) 


In particular for the interaction potential of “semi-transparent sphere” type (6.14) , we have: 

1 

U 3 ,(r,r') = U 3m 0(r 3 - r)0(r 3 -r')(2Z + l)[0(r 3 - r-r')d/ o + 0(r+ r'-r 3 )^ JPi(x) dx], (8.3) 


P(f) = P(°) = f = ^ 2 - pifa,r') = ^^jdkk 2 ji(kr)ji(kr’). (8.4) 


/c 3 1 

ai(r) = U 3m d(r 3 -r )—^~— 
2 ^ J fc F r 


« 2 (r) = ^^(fc F r 3 ) 3 0 (r 3 - r) 
7Zjt 


k^rs 


j$VcFr 3 )-j 0 (k ¥ r 3 )j 0 [k ¥ (r 3 - r)] + J dyj 0 (A: F r-y);i(y) 


— -12 — +16 


-fc F (r 3 -r) 


f3 


+3 


5^-2 5 

b\ = — p 2 ( 0 )[/ 3 m r 3 6 =—^t/amtfcpra ) 6 = 0.0023(/ 3m (fc F r 3 ) 6 , 


^2 = 


U 3lt 


B 3 {k^r 3 ) - 12 (Ar F r 3 ) 5i(fc F r 3 ) + 16(fc F r 3 ) -Bo(fc F r 3 ) , 


12tt 2 

where D ra (z) = Jq y n jf (y)dy. When calculating «i(r), the summation formula is used (see 0D> P-133) 


(8.5) 


00 I /- 1 

£(21 + l) 7 /(u) 7 ;(n)P;(x) =jo V it 2 + v 2 -2uvx . 

1=0 v ' 


(8.6) 


In the limit fc F r 3 « 1, we have more simple formulae 


b 2 ‘ 


216n 2 


U 3m {k F r 3 f, «i(r)i 


U 3n 


432tt 3 7 


r (fc F r 3 ) fa 0 (r 3 - r) 


+3 


— -12 — +16 


r 3 


(8.7) 


Note that even at fc F r 3 = 1, the formulae ( 8 .7) give a mistake of the order and less than 10%. For the follow¬ 
ing, it is convenient to define a characteristic value of density through the radius of the pair interaction: 


1 _ 47T 3 

rc* 3 2 


( 8 . 8 ) 
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If, for example, r 2 = 3 x 10 -8 cm is taken, then n t - 0.88x 10 22 cm -3 . This density is close to the particle 
number density in liquid. Note that, for example, the particle number density in liquid helium-3 at near¬ 
zero temperature and vapor pressure: m He ~ 1.6xl0 22 cm -3 . Since fcpr 3 = (97r/4) 1/3 (r 3 /r 2 )(rc/n*) 1/3 , then 
the condition fcpr 3 « 1 can be written in the form 

(r 3 / r 2 ) («/«*) 1/3 « 1. (8.9) 


Since it should be considered likely that r 3 ~ r 2 , the noted condition is satisfied in the limit of low den¬ 
sities. However, as was noted, even at n ~ n, the formulae {8.7} are valid with good accuracy. Let us 
introduce the designation that will be used below for the dimensionless density 



( 8 . 10 ) 


where r s = [3/(4 tt7z)] 1/3 is the radius of a sphere whose volume equals the volume per particle. It is im¬ 
portant to note that, even in the case when the dimensionless density (8.10} is considerably less than 
unity, the value of the dimensional density proves to be much greater than that of low-density gases. In 
the following formulae we will also use the next designations: 


r 3 

do 

U 3 m 

- y 

a = —, 

U = - 

r 2 

r 2 

u 2m 


( 8 . 11 ) 


<%o = m(7 2m r|/3lJ 2 . With account of 1 7.9', the relations 8.4 , d8.s| determine the contribution of three 


body interactions into the direct and exchange parts of the self-consistent potential 


W 0 (3) = 2 foi - b 2 , W' 3) (r) = ai(r)-3/i^^ 

u L fcpr 

The contribution of pair forces into these quantities for the potential of 
(6.12} is given by the formulas 


a 2 {r). (8.12) 

“semi-transparent sphere” type 


= U 2m n. 






A (for) 

kpr 


0 , 


r < r 2 , 
r > r 2 . 


(8.13) 


The obtained relations (8.12} , (8.13} for the self-consistent potential enable us to determine by the 
formula (7.12} the contribution of pair and three-body interactions into the effective mass. With account 
of only pair forces, the effective mass is determined by the expression 


m 6 a 

- — 1 H- - -~ 

ra * 2 n {k F r 2 r 


B 2 {k^r 2 ). 


(8.14) 


Hence, it follows that for the repulsive pair interaction U 2m > 0. the effective mass of a quasiparticle 
proves to be less than the mass of a free particle, and for the attraction U 2m < 0, it is greater. The effective 
mass with account of three-body interactions in this case has the form 


m 
m * 


m 

m* 2 


r 3 

4nm r o \2nmn C 2 2 

^2 J ai[r)j\{k P r)r dr + J j\ (k F r)a 2 (r)r dr. 

o F o 


(8.15) 


It should be stressed that, for all densities at which the self-consistent field approximation holds, the 
formula (8.15} is exact. At low relative densities, when the condition (8.9} is satisfied, the formula for the 
total effective mass (8.15} can be represented in the form of expansion in powers of density 


m 3 3 I9tt V 

- = In- an -— 

m* 10 701 4 j 


2/3 


-5/3 18 -2 

a n h- auX n 

2560 


= 1 + 0.3 a n - 0.158 a n 513 + 0.0621 a uX 3 n 2 . 


(8.16) 
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Figure 1. (Color online) Dependencies of the effective mass on density: (1) a = -2.5, u = 0; (2) a = -2.5, 
u = -0.3; (3) a = 2.5, u = 0; (4) a = 2.5, u = 1. It is everywhere assumed that A = 1. 


In (8.14) , (8.16) and below we use the designations (8.11) . The second and the third terms in (8.16) are 
determined by pair forces, and the fourth by three-body forces. As it is seen, with increasing density the 
contribution of three-body forces into the effective mass rises faster than that of pair forces, and the 
relative role of three-body forces rises. Furthermore, attention is drawn to a strong dependence (as A 8 ) 
of the contribution of three-body forces on the ratio of the radii of three-body and pair interactions. 

Dependencies of the effective mass on density, calculated by the exact formula (8.15) , are shown in 
figure [l] For the attractive interaction and with neglect of three-body forces, the effective mass proves 
to be greater than the mass of a free particle and monotonously rises with an increasing density, reach¬ 
ing its maximum. With a further increase of density, the rise of the effective mass changes into its fall, 
although, as before for all physically reasonable densities, it remains greater than the mass of a free 
particle (curve 1 in figure]!]. Accounting for the three-body interaction with a positive constant leads to a 
reduction of the region of the rise of the effective mass and to its faster decrease at high densities (curve 2 
in figure]!}. For the repulsive interaction and in the absence of three-body forces, the effective mass for 
all reasonable densities is less than the mass of a free particle (curve 3 in figure]!}. Accounting for the 
three-body interaction with a positive constant leads to a faster monotonous decrease of the effective 
mass (curve 4 in figure]!}. _ 

The total pressure p - -D./V of the Fermi system with the three-body interaction of the form 1 6.10 
at finite temperatures, according to the general formulae (5.3), (5.4), is as follows 


P 

Po 


P3 


P0 + P2 + P3, 


2 T 

-^,2 


h 2 k 2 \ 

P: 


2 nu 


P 2 = 4 n 


UU <_xj 

- J U 2 {r)p 2 [r) r 2 dr + 2p 2 (0) J U 2 {r)r 2 dr 


J drdr' U 2 {r, r', |r- r'|)[- 4p(0)p 2 (r') + ^p 3 (0) + ^p(r)p(r')p(|r- r'|)], (8.17) 


where the term p 0 is a contribution of a gas of fermions with the effective mass m„ and the terms p 2 and 
P 3 give, respectively, contributions of the pair and three-body interactions. Together with the formula 
(7.17) for the particle number density, formulae (8.17) define in a general form the system’s equation of 
state for the considered form (6.10) of three-body interactions. 

At zero temperature, the pressure of the fermion gas 


Po — 


(3tt 2 ) 2/3 h 2 


5/3 


(8.18) 


and contributions of the pair and three-body forces into the total pressure (8.17) in the model of “semi¬ 
transparent sphere” potentials for them [formulae (6.12), (6.14)], with the expansion formulae (8.1)-(8.4) 
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taken into account, acquire the form 


P 2 = 


Uim n * 


[- nB 0 {k F r 2 )- n 2 \, 

L 71 j 


CXJ 

(2 \ 167T r o 

p 3 = 4p(0)^-i>i-h 2 J + — I p(r)ai(r)r 2 dr, 


(8.19) 

( 8 . 20 ) 


where b\, b 2 , a\ (r) are defined by the formulae (8.5). In the case of densities using the expansion 


7 72+3 


-,n +5 


B n {z) i 


9(rc + 3) 45(n + 5) 525(ti + 7) 

we find the dependence of the pressure due to the pair interaction on density 
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( 8 . 21 ) 


Here the pressure is related to that of a gas of particles at the characteristic density |8.8) : p 0 * = 
\(3n 2 ) 213 [ti 2 1 m)n s J 3 . The sign of pressure 18.211 is determined by the sign of the pair interaction con¬ 
stant and, at negative value of this constant, the interaction contribution into the pressure is negative. 
The main term in the expansion of pressure due to the three-body interaction (8.20) has the form 
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a uA 8 n 11/3 = 0.186a u A 8 fi 11/3 . 


( 8 . 22 ) 


Attention should be paid to the fact that the expansion of p 2 holds in even powers of the quantity n 113 
and the expansion of p 3 in odd powers, and the latter begins with a high power. 

Since the effective mass depends on density, then the pressure of a gas of quasiparticles (8.18) , with 
account of (8.16), can also be represented in the form of expansion in powers of density: 
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2560 


= n 513 + 0.3 a h 813 - 0.158 a fi 10/3 + 0.0621 a wA 8 fi 11/3 . (8.23) 

Taking into account (8.21)-(823), we find the expansion of the total pressure in powers of density: 
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(8.24) 


Comparison of the accurate dependence of the total pressure on density, calculated by the formulae 
(8.18) (8~20) and the approximate dependence, calculated by the formula (8.24) , is given in figure [2] It 
is seen that even at n ~ 1, the pressure, calculated by the approximate formula, differs weakly from its 
accurate value. 

Let us discuss the issue of the thermodynamic stability of the Fermi system, assuming h < 1. For a sys¬ 
tem to be stable, it is necessary that its compressibility or (which is the same) the squared speed of sound 
should be positive, so that dpldn> 0. Since in this case the two first terms give the main contribution into 
the total pressure (8.24) and qualitatively taking into account the new effect of three-body interactions, 
we find the condition of stability in the form 

9 ( 4 ) 2/3 m 1749 „ o 

1+- — afr 3j t -az^A 8 fr>0. (8.25) 
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Figure 2. (Color online) Comparison of the accurate dependence of the pressure (1), calculated by the 
formulae |8.18) - (Oo) and the approximate dependence (2), calculated by the formula (8.24) , at a = -0.3, 
u = -1, X - 1. (3) Ideal Fermi gas. 


For positive constants of pair and three-body interactions, the system is always stable. More interesting 
is the case when the constant of the pair interaction is negative. Then, without account of three-body 
interactions, the condition of stability should be satisfied which can be represented in equivalent forms 


n m < 


1/3 


lal -1 ~ 0.82|ar 1 , 


/18 ) 1/3 

/"s > I ^ I |a 0 l ~ 1.22 |a 0 |. 


(8.26) 


Accounting for the three-body interaction, if it is of a repulsion character, extends the region of stability 
and can lead to stabilization of the system with the negative pair potential at arbitrary densities in case 
of fulfilment of the following (8.25) condition 


\u\ > 


200 5 5 \a\ 5 
583 7T 4 I 8 


0.408 
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(8.27) 


Some dependencies of the pressure on density for different signs of the interparticle interaction are 
shown in figure [3] In the case of attraction and with neglect of three-body interactions, the spatially 



n 

Figure 3. (Color online) Dependencies of the pressure on density: (1) a = -0.3, u = 0; (2) a = -0.3, u = 
-0.5; (3) an ideal Fermi gas; (4) a = 0.3, u = 1. It is everywhere assumed that A = 1. 
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Figure 4. (Color online) The effect of three-body repulsive forces (f/ 3 m > 0) on stability of the Fermi 
system with pair attractive forces: (1) u = 0; (2) u = -0.085; (3) u = -0.098; (4) u = -0.11; (5) u = -0.13. It 
is everywhere assumed that a = -0.3 and X = 1. 


uniform state is stable only at low densities for which the condition (8.26) is satisfied. At high densities, 
the pressure decreases with an increase of density (curve 1 in figure [3j and the spatially uniform state 
ceases to be stable. Sufficiently intensive three-body repulsive forces lead to stabilization of the system 
with pair attractive forces (curve 2 in figure[3}. 

The effect of three-body repulsive forces on the stability of the Fermi system is illustrated more in 
detail in figure [4] With an increasing strength of the three-body repulsive interaction, the regions of 
stability of such a system extend (curve 2, 3 in figure [4), and a further growth of the intensity of the 
three-body interaction leads to stabilization of the system for all physically reasonable values of density. 

The form of the quasiparticle energy spectrum with account of the interparticle interactions is shown 
in figure [s] Accounting for only the pair attraction leads to a slowing down of the energy increase with 
increasing momentum relative to an ideal Fermi gas (curve 2 in figure [5). Accounting for the repulsive 
three-body forces, against the background of the attractive pair forces, additionally yields a small increase 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

K 

Figure 5. (Color online) The quasiparticle energy spectrum e(jc) (in the variables e = e^./ep, k = fc/fcp): 
(1) Ideal Fermi gas [e(k) = k 2 - 1, m* = m*lm = 1]; (2) a = -2.5, u = 0, n = 1 (m* = 1.78); (3) a = -2.5, 
u = —1, X - 1, h = 1 (w* = 1.50); (4) a = 2.5, u = -1, X = 1, n = 1 (m* = 0.75). 
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of the quasiparticle energy (curve 3 in figure[s}. The effect of three-body forces on the spectrum is appre¬ 
ciably weaker than the role of pair forces, so that accounting for the attractive three-body forces in the 
presence of the pair repulsion weakly prevents the increase of the quasiparticle energy as momentum 
increases (curve 4 in figure[s}. 


9. Conclusion 

The self-consistent field equations are obtained, and, within this model, thermodynamic relations are 
derived for a normal system of Fermi particles with account of both pair and three-body forces. The 
satisfaction of the essential requirement of fulfilment of all thermodynamic relations already in the self- 
consistent approximation leads to a unique formulation of the self-consistent field model. Emphasized 
are the novelty and universality of the proposed approach to the formulation of this model and the use¬ 
fulness of the method for describing not only Fermi, but also Bose systems, in particular phonons, and 
relativistic quantized fields. The approach being developed allows us also, as shown in the paper, to ac¬ 
count for many-body interactions naturally. 

It is shown that three-body interactions of zero radius give no contribution into the self-consistent 
field and, in order to account for the effects due to such interactions, it is necessary to account for their 
nonlocality. The case of the spatially uniform system is considered in detail. General formulae are derived 
for the system’s equation of state and the effective mass of quasiparticles with account of three-body 
forces. Dependencies of the quasiparticle effective mass and the system’s pressure on density at zero 
temperature are obtained, with pair and three-body forces accounted for in the model of interaction 
potentials of “semi-transparent sphere” type. 

It is shown that pair interactions of repulsive character reduce the quasiparticle effective mass rela¬ 
tive to the mass of a free particle while attractive pair interactions, on the contrary, raise it. The effective 
mass and pressure are numerically calculated at zero temperature and expansions of these quantities are 
derived in powers of the relative density with account of three-body forces. It is shown that the relative 
contribution of three-body interactions into thermodynamic quantities rises with an increasing density. 
The effect of three-body forces on the stability of the Fermi system is considered, and it is shown that in 
the case of repulsion, their being taken into account extends the region of stability and can lead to stabi¬ 
lization of the system with pair attraction. The quasiparticle energy spectrum is calculated with account 
of the interparticle interactions. 
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Mo^ejib caMoy3rofl>i<eHoro nosin j \ j\r (j>epivii-CMCTeM 
3 BpaxyBaHHflM Tpi/macTi/iHKOBi/ix B3aeM0flm 

K).M. rio/iyeKTOB, O.O. CopoKa, C.M. LLIy/ibra 

iHCTHTyT TeopeTHHHoi <t>i3MKH iM. O.I. Axie3epa, HHLf XOTI, By/i. AxafleMiHHa, 1, 61108 XapKiB, yKpa'ma 

Ha ocHOBi MiKpocKoniHHoi MO.ge.ni caMoy3roflxeHoro no/ia no6yflOBaHO TepMOflHHaMiKy cncreMH 6araTt,ox 
tjjepMi-MacTMHOK npn CKiHHeHHX TeMnepaTypax 3 BpaxyBam-iflM TpHHaCTHHKOBHX B3aeMO,giH i OTpuMaHO piBHfl- 
HHfl pyxy KBa3i4acTHH0K. noKa3aHO, mo fle/ibTanofli6Ha TpHnacTHHKOBa B3aeMO,gm He flae BHecxy b caMoy3ro- 
flxeHe no/ie, i pp in onncy TpHHaCTHHKOBHX ch/i cnifl BpaxoByBaTH ix He/ioKa/ibHiCTb. /feTa/ibHO po3r/iBHyTO npo- 
CTopoBOOflHopiflHy CHCTeMy i b paMKax po3BHHeHoro MiKpocKoniHHOro niflxofly OTpuMaHO 3ara/ibHi <|)opMy/in 
fl/lfl ec|>eKTHBHOi MaCH 4>epMiOHy i piBHflHHfl CTaHy CHCTeMH 3 BpaxyBaHHBM BHeCKy TpHHaCTHHKOBHX B3a£MOfliii. 
/f/iB noTeHpiaay Tuny "HaniBnpo3opoT cc()epn" npn Hyyii TeMnepaTyp nnce/ibHO po3paxoBaHO e<j>eKTHBHy Macy i 
tuck. 3HafifleHO po3BHHeHHB ecfieKTnBHOi Macn i Tncxy 3a cryneHHMH Lpi/ibHOCTi. riOKa3aHO, m° npn BpaxyBaHHi 
Ti/ibKn napHHX ci/iti, B3aeMOflin : m° Mae xapaKTep BiflLUTOBxyBaHHn, 3MeHLuye etj)eKTHBHy Macy KBa3iHacTHHKH 
nopiBHBHO 3 Macow Bi/ibHoT HacTnHKn, a y pa3i npnmraHHB — 36i/ibiuye. Po3mnHyTO nnTaHHn TepMOflHHaMi- 
hhoT CTiiiKOCTi ((jepMi-cncTeMn i noKa3aHO, mo TpnnacTHHKOBa B3aeMOflia, axa Mae xapaKTep BiflLUTOBxyBaHHfl, 
p03LUHpK)e 06/iaCTb CTiHKOCTi CHCTeMH 3 MiXCHaCTHHKOBnM napHHM npHTBTaHHBM. Po3paXOBaHO eHepreTHHHHH 
CneKTp KBa3inaCTHHKH 3 BpaxyBaHHBM TpHHaCTHHKOBHX CH/1. 

IOuohobI c/ioBa: caMoy3rop,>KeHe none, Tpi/inacTHHKOBi B3acM0flii, ecpeKTMBHa Maca, cpepiviiOH, 
piBHBHHB CTaHy 
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